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ABSTRACT 

We simulate Gamma-Ray Bursts arising from internal shocks in relativistic winds, calculate their power density 
spectrum (PDS), and identify the factors to which the PDS is most sensitive: the wind ejection features, which 
determine the wind dynamics and its optical thickness, and the energy release parameters, which give the pulse 
50-300 keV radiative efficiency. For certain combinations of ejection features and wind parameters the resulting 
PDS exhibits the features observed in real bursts. We found that the upper limit on the efficiency of conversion 
of wind kinetic energy into 50-300 keV photons is ~ 1%. Winds with a modulated Lorentz factor distribution of 
the ejecta yield PDSs in accord with current observations and have efficiencies closer to 10~ 3 , while winds with 
a random, uniform Lorentz factor ejection must be optically thick to the short duration pulses to produce correct 
PDSs, and have an overall efficiency around 10~ 4 . 

Subject headings: gamma-rays: bursts - methods: numerical - radiation mechanisms: nonthermal 



1. INTRODUCTION 

Internal shocks occurring in a transient, unstable relativistic 
wind (Rees & Meszaros 1994) are believed to be the source 
of the ^100 keV emission and the complex temporal struc- 
ture of Gamma-Ray Bursts (GRBs). The shocks resulting from 
such instabilities heat the expanding ejecta, amplify pre-existing 
magnetic fields or generate a turbulent one, and accelerate elec- 
trons, leading to synchrotron emission and inverse Compton 
scatterings. Synchrotron self absorption and pair formation 
may further shape the emergent spectrum, depending on the 
choice of model parameters, as shown by Papathanassiou & 
Meszaros (1996) and Pilla & Loeb (1998). The spectrum 
of the emitted radiation has also been analyzed by Daigne & 
Mochkovitch (1998) and Panaitescu & Meszaros (1999). The 
synchrotron self-Comptonized emission from internal shock 
GRBs has also been studied by Sari & Piran (1997) and Ghis- 
ellini&Celotti(1999). 

In this paper we analyze the temporal features of GRBs 
obtained in the framework of relativistic winds, through the 
burst power density spectrum (PDS). The recent work of Be- 
loborodov, Stern & Svensson (1998) has put into evidence in- 
teresting PDS features, which we use for analyzing some of the 
model parameters and the features of the wind ejection. We use 
the observed integral burst intensity (log N — log P) distribution 
(Pendleton et al. 1996) as a constraint in our choice for some 
of the model parameters. Other properties of the observed light 
curves, such as the distribution of pulse durations (Norris et 
al. 1996) or those of pulse fluences and time intervals between 
peaks (Li & Fenimore 1996), can also be used in the study of the 
wind ejection features (Spada, Panaitescu & Meszaros 1999). 

2. OUTLINE OF THE MODEL FEATURES 

The model used here is similar to that developed by Daigne 
& Mochkovitch (1998), and is different in the following ways. 
Our treatment of the radiation emission takes into account the 
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up-scattering of the synchrotron photons, which may be the 
dominant emission process at ~ 100 keV if the magnetic field is 
sufficiently low. The optical thickness (to Thomson scatterings) 
of the wind is taken into account, which may be very important 
for the burst PDS and overall efficiency. We use the shock jump 
equations to determine the physical conditions in the shocked 
fluid and we include the effects arising from electron cooling in 
the shape of the emission spectrum. 

After setting the dynamics of the wind ejection, we calcu- 
late the radii where internal collisions take place, determine the 
relevant physical parameters in the shock fluid - bulk Lorentz 
factor (LF), typical electron random, magnetic field, etc. - and 
calculate the features of the emitted radiation: energy of the syn- 
chrotron and inverse Compton spectral peaks (with allowance 
for the Klein-Nishina effect), Comptonization parameter, 50- 
300 keV emission (accounting for scattering during propagation 
through the wind), and its duration. These quantities are neces- 
sary for the computation of the observer frame pulse features: 
fluence, duration, arrival time, taking into account relativistic 
and cosmological effects. Assuming the two-sided exponential 
pulse shape identified by Norris et al. (1996) in real bursts, one 
can calculate the peak photon flux for each pulse. Below we 
describe the most important aspects related to the calculation of 
the some of the above mentioned quantities. 

The wind is discretized as a sequence of shells, released by 
the central source during a time t w and with an average interval 
t v between consecutive ejections, so the number of shells is just 
N = t w /t v . The mass of each shell is drawn from a log-normal 
distribution determined by the average shell mass M — M w /N, 
M w being the wind mass, and a dispersion <tm = M. The LF 
r\ of each shell is randomly drawn from the interval [rj m , t/m]- 
Both T] m and t)m can be chosen constant during the entire wind 
ejection, so that r\ has a uniform distribution (we shall refer 
to this case as the "uniform wind"), or can vary on time-scale 
comparable to t w and much larger than t v , according to a certain 
law. For brevity, in the latter case, we shall keep rj m constant 
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and use an tjm that varies with ejection time as a square sine 
with only one oscillation (this will be the "modulated wind"), 
thus the LF of the j lh ejected shell is 

Vj =Vm + aj sin 2 (irj/N)(i] M - Vm) 1 <3 < N , (1) 

where a, ; is a random number between and 1 . The time elapsed 
before the ejection of a shell is chosen to be proportional to the 
energy of that shell, thus leading to a wind luminosity L w that 
is constant on average throughout the entire wind ejection. 

The shock-accelerated electrons have a power-law distribu- 
tion of index —p, starting from a low random LF set by the 
electron injection fraction ( < 1 and the fraction e e of the 
internal energy stored in the electrons: 
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where v! and n' e are the comoving frame internal energy density 
and electron number density of the shocked fluid, respectively, 
which we calculate with the aid of the shock jump conditions, 
and Cdis is the dissipation efficiency, i.e. the fraction of the 
kinetic energy that is converted into internal. The magnetic field 
is assumed turbulent and parameterized through the fraction Eb 
of the internal energy density of the shocked gas stored by it: 
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For the computation of the burst emission in the 50-300 ke V 
band, we approximate the synchrotron spectrum of each pulse 
as a three power-law function, with slopes that depend on p 
and on the relative values of the cooling and peak frequencies. 
A fraction min(l,T£) of the synchrotron photons undergoes 
max(l , t 2 ) up-scatterings in which their energy is increased by 
a factor ~ (4/3)7 2 „ per scattering (unless the Klein-Nishina 
regime is reached), where tq = <JTh{C, n 'e) nun (ct', A') is the 
optical thickness to hot electrons, A' being the co-moving shell 
thickness, ti ( — t' sy /{\ + y) the radiative cooling time-scale, 
t'sy x {lm,B 2 )~ l the synchrotron cooling time-scale and y the 
Comptonization parameter, which depends on (for < 1, 
V = lm T 0- The A' is calculated from the shock compression 
factor and from the thickness of the shell before the collision. 
For the latter we assume that during the free expansion between 
successive collisions, the shell thickness fractional increase is 
equal to the fractional increase of its radius r (d log A'/d log r = 
1). 

We take into account the optical thickness to Thomson scat- 
terings on cold electrons in the emitting shell (tq) and in the 
outer part of the wind through which a pulse propagates (r), 
by diminishing the pulse intensity by a factor (1 — e -r °)TQ 
and e _T , respectively. However we do not include in the burst 
light-curve calculation the scattered photons which, in the end, 
still arrive at the observer. If the cumulative time delay due 
to scattering/diffusion through the wind is larger than the burst 
duration then, for the observer, the photon scattering practi- 
cally acts like an absorption process. Otherwise it may mimic 
an emission process with a timescale shorter or comparable to 
the burst duration, and one should take it into account in the 
light-curve (and PDS) calculation. 

Very important for the burst PDS is the computation of the 
pulse duration, which is determined by: 



1) the geometrical curvature of the emitting shell (the observer 
receives most of the radiation from a zone extending up to an 
angle ~ T~ 1 relative to the central line of sight), 

2) the lab-frame electron radiative cooling time t 1 = Ttiy, 

3) the lab-frame shock's shell-crossing time t& — A/|u s ^ — vq\, 
v s h being the shock's speed, determined from the hydrodynam- 
ics of the collision, and vo the shell pre-shock flow velocity. 
The spreads in the photon arrival time are functions of the angle 
relative to the central line of sight; so is the intensity of the rel- 
ativistically boosted emission. The intensity-averaged observer 
time spreads corresponding to the three factors above are 
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which we add in quadrature to calculate the pulse duration 
ST. For further qualitative estimations we will consider that 
the pulse duration is positively correlated with the radius where 
the collision takes place. This would be obvious if ST ~ STg, 
but it is a correct assumption because both <5T 7 and 5T& have 
the trend of increasing with r, due to the fact that the elec- 
tron LFs/magnetic fields decrease with radius, while the shell 
thickness increases. 

The 30-500 keV pulse fluence is a fraction of the kinetic 
energy of the colliding shells, equal to the product of 

1) the dissipation efficiency edis- If Vm 3> Vm then edis can 
exceed 10-20% for the collisions taking place at small radii, 
where there is a larger difference between the shell LFs, de- 
creasing to 1% or less for the late collisions. 

2) the radiative efficiency at which the internal energy is con- 
verted into radiation. It is upper bounded by e e , the fractional 
energy in electrons, reached if the 7 m -electrons are radiative. 

3) the window efficiency £2+3, representing the fraction of the 
radiated energy that arrives at observer in the 50-300 ke V band 
(the 2 nd and 3 rd BATSE channels). Given the broadness of the 
synchrotron and inverse Compton spectra, £2+3 cannot exceed 
30%forp = 2.5. 

Thus the overall burst efficiency can hardly exceed 1 %, a value 
reached if the model parameters are such that for the most en- 
ergetic collisions in the wind the above efficiencies are close to 
their maximal values. 

3. NUMERICAL RESULTS 

Beloborodov et al. (1998) (hereafter BSS98) have calculated 
the average PDS of 214 GRBs longer than 20 seconds using 
their 50-300 keV light-curves with 64 ms resolution. The PDS 
they obtained follows a power-law Pt oc /~ 5 / 3 (/ is frequency) 
from /; ~ 0.03 Hz to fh ~ 1 Hz, decaying weaker below /; and 
stronger above In other words, / 5 / 3 P/ is flat up to fh, where 
a break is observed. We have identified three factors that can 
lead to such PDSs: i) the pulse window efficiency £2+3, H) the 
LF ejection law, which determines the wind dynamics and the 
dependence of e<& s on ST, and Hi) the wind optical thickness r 
for short pulses. 

Figure la shows the effect of C an d Eb on e pu i se , the frac- 
tion of the wind entire kinetic energy that is radiated by each 
pulse in the 50-300 keV band. The edis and n' e decrease with 
r, thus 7,„ (eq. [^]) and B (eq. [0]) exhibit the same trend, 
which makes £2+3 decrease with ST. The two energy release 
parameters ( and e b alter the peak energies of the synchrotron 
and inverse Compton spectra (they are higher for lower Q and 
the Comptonization parameter (which is larger for smaller eb)- 
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In the case shown here, decreasing these parameters results in 
increasing the £2+3 of the longer pulses. The PDSs in Figure 
16 illustrate how the burst power is distributed versus variabil- 
ity timescale, by combining the pulse duration distribution with 
the e pu i se -ST dependence of Figure la. The conclusion to be 
drawn from Figure 16 is that, for uniform winds that are op- 
tically thin to most pulses, electron fractions below unity and 
magnetic fields well below equipartition reduce the power in the 
short pulses, thus representing one possible cause for the 1 Hz 
break; nevertheless the PDS remains flat below 1 Hz. 

An ejection with modulated LFs mitigates the decrease of 
with radius, by clumping shells earlier in the wind expansion and 
producing groups of shocked shells that travel longer distances 
before suffering more collisions. Thus the differences between 
the LFs of the shells colliding at larger radii are greater, yielding 
a better e pu i se for longer pulses. This is illustrated in Figures 
lc and Id for winds that are optically thin, where we used a LF 
ejection modulated by a sine square (eq. [jj]]). Note the flatness 
of / 5/3 P/ at / < 0.5 Hz. 

The fact that ST and the radius where the collision takes place 
are positively correlated suggests that the observed lack of power 
in short duration pulses may also be due to the large r of the wind 
at early lab-frame times. We shall refer to the case where short 
pulses occur predominantly below the photospheric radius as the 
"optically thick case", without meaning that the wind is thick for 
all pulses. Equation (3) from Rees & Meszaros (1994) gives an 
estimation of the photospheric radius r p h — (nL w )/(87rc 3 77 3 ), 
where k = 0.4cm 2 g~' and rj is the wind average LF. If ST is 
mainly determined by the geometrical curvature of the emitting 
shell, then the duration of the pulses emitted at the photospheric 
radius is ST ph ~ r ph /(2rj 2 c) : 

kL 

ST ph = 77 — 7=7 = l-QLwMVi 5 s • ( 5 ) 
Mmcrr) 

Equation (|J) shows that the duration of the pulses for which 
the wind is optically thick is very sensitive to rj. An increase in 
the wind's r may result from a higher L w or a lower rj, either 
case implying a higher wind mass. Numerically we found that 
for an uniform wind with rj m = 30 and L w = 10 54 erg s _1 t\m 
must be below 300 to increase the wind's r to the point where 
there is substantial loss of power at the high frequency end of 
the PDS. Figures le and If show that this loss of power can be 
enhanced if the magnetic field is far from equipartition. 

To account for the effects arising from the cosmological dis- 
tribution of GRBs in the calculation of the PDS and intensity 
distribution, the GRB redshifts are chosen from the probability 
distribution 

dp n c (z) dV 

7*7- — 77-) (0) 

dz 1 + z dz 

where dV/ dz is the cosmological comoving volume per unit 
redshift, and n c (z) is the rate density evolution of GRBs. For 
brevity, we consider here the case of a constant rate density, 
keeping in mind that other reasonable functions n c (z) alter the 
PDS, though the changes are not drastic. We use an un-evolving 
power-law distribution for the wind luminosity: 

3>(L W ) oc L-P , L m < L w < L M , (7) 

and zero otherwise. Note that this not the same as assuming 
that the GRB 50-300 keV luminosity at the source redshift has 
a power-law distribution, as it is usually done (e.g. Reichart & 
Meszaros 1997, Krumholz et al. 1998). 



Figure 2 illustrates the effect of L w and z on the PDS, for 
a sin 2 -modulated wind. As can be seen, there is a shift of 
power from higher frequencies to lower ones if L w is increased. 
The same is true for a burst placed at a higher redshift. The 
latter is due to the time dilation, which makes pulses appear 
longer. An increase in L w leaves unaltered the wind dynamics 
and enhances the magnetic field (the comoving particle density 
in eq. [||] is higher). Thus the emission becomes harder and 
the window efficiency will favor longer pulses. If the wind is 
quasi-optically thick, the correlation between the PDS and L w 
is strengthened by the increase of r with L w . 

Figure 3a shows the intensity distribution of simulated bursts, 
using the peak fluxes on 256 ms timescale, compared to the ob- 
served distribution (Pendleton et al. 1996). About 50% of the 
400 simulated bursts are brighter than 0.4 7cm _2 s~ 1 , thus the 
number of bursts in our sample is close to that of the bursts an- 
alyzed by Pendleton (1996). Excluding the two points with the 
lowest peak flux and the one with the highest (which contains 
only 1-3 bursts), the \ 2 f° r the sets of optically thin, modulated 
and thick, uniform winds are 7.8 and 16.4, respectively, for 8 
degrees of freedom. Thus only the intensity distribution of the 
modulated wind bursts is consistent with the observations. A 
lower L m would yield a better fit of the log N — log P distri- 
bution for uniform winds, but the PDS would have too much 
power at high /. 

The PDSs shown in Figure 36 have been calculated by av- 
eraging the PDS of all bursts with peak photon fluxes brighter 
than 1 7 cm~ 2 s~ 1 on the 64 ms timescale. More than one of the 
mechanisms presented in Figure 1 must be invoked in order to 
obtain /~ 5 ' 3 behavior of the PDS and the 1 Hz break: for mod- 
ulated bursts we used £ = 10~ 2 and eb = 10~ 6 ; for uniform 
winds ijm = 220 and eb — 10 -4 reduce the brightness of the 
short pulses to the point where f 5 ^ 3 Pf exhibits a flat part, albeit 
over a range of frequencies narrower than observed. In the for- 
mer case the wind efficiency is slightly larger than 0. 1 %; in the 
latter, the scatterings of the photons in the wind diminishes the 
burst efficiency to ~ 10~ 4 . For this reason, thick un-modulated 
winds require a larger energy budget than the thin modulated 
ones. 

Figure 3c shows how the power spectra of individual bursts 
are spread around the average PDS. The distribution was cal- 
culated by averaging the Pf distributions obtained for a large 
number of frequencies. This figure also shows the exponential 
fit dN/dlogP f = (In 10)(Pf/P f )exp(-P f /P f ) that BSS98 
found to approximate well the observed distribution. Note that 
the two models for ejection LF partly bracket the exponential fit, 
suggesting that in real bursts the wind is neither as modulated 
as the squared sine employed here, nor as erratic as the uniform 
distribution, but somewhere in between. This conclusion is con- 
sistent with the fact that the observed 1 Hz break is in between 
the breaks exhibited by the simulated PDSs shown in Figure 36. 

In Figure 4 we show the typical 50-300 keV light-curves of 
bursts with a modulated and an uniform wind. 



4. CONCLUSIONS 

We have investigated the power spectra of the 7-ray light 
curves expected from internal shock GRBs. The PDS is found 
to be most sensitive to the wind luminosity L w , its average LF 
v(Vm, Vm), an d the burst redshift z. The PDS is sensitive to 
other model parameters, such as the wind duration t w , if the 
wind LF is modulated, and the variability timescale t v , if the LF 
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distribution is just random. The properties of the PDS are also 
affected by the dynamics of the unsteady relativistic wind and 
by the microscopic parameters of the shocked fluid. 

We have compared our results to those of Beloborodov et 
al. (1998), and have identified three possible reasons for the 
observed deficit of short (ST < 1 s) pulses: 

1) the 50-300 keV radiative efficiency of such pulses may be 
smaller than for longer pulses, due to low electron injection 
fractions and magnetic fields well below equipartition. 

2) the short pulses may result from the collision of light shells, 
carrying little energy, due to a modulation of the ejection LFs. 
Winds with modulated ejection LFs (here we used a squared 
sine) yield PDSs and log N — log P distributions consistent with 
the observations, and have a 50-300 keV efficiency <; 10~ 3 . 

3) the short pulses may occur below the photospheric radius. 
The duration of pulses occurring at the photospheric radius 
depends strongly on the average LF of the wind (eq. [|[]). For 
a wind LF distribution that has a lower bound r\ m = 30, an 
upper limit r)M — 220 leads to a significant decrease of the 
burst power at high frequency, yielding a PDS behaving like 
/~ 5 / 3 for 0.2 Hz < / < 2 Hz, and to a 50-300 keV efficiency 
around 10~ 4 . A lower t)m yields a PDS with even lower power 
at high frequency, but the wind efficiency in the middle BATSE 
channels becomes too small. 

The overall burst efficiency is small, due to the low dissipation 



efficiency of the wind, typically between 1% and 10%, and 
to the broadness of the synchrotron self-Comptonized spectra 
from power-law distributions of electrons, leading to window 
efficiencies ~ 10%. For the wind luminosities and durations 
used in the calculations whose results are shown in Figure 3, a 
beaming factor ranging from ~ 10 (for bursts with L w = L m ) 
to ^ 1000 (for L w — Lm) is needed to maintain the energy 
requirements below the upper limits found by Meszaros , Rees & 
Wijers ( 1 999) for the energy that can be extracted from plausible 
GRB progenitors. 

The study of the power spectra of internal shock GRBs has 
shown that if the ejection parameters of optically thin winds are 
totally random, the resulting spectrum is flat, with equal power 
at low and high frequency. In order to explain the observed 
/~ 5 / 3 behavior of the PDS, the wind must be modulated such 
that collisions at large radii release more energy in the observing 
band than the collisions occuring as small radii. A modulation 
of the wind ejection is physically quite plausible, and the fact 
that it is necessary to introduce this in order to obtain the correct 
PDS is one indication of the value of the power spectrum as a 
tool in studying the physics of the GRB "central engine". 

This research is supported by NASA NAG5-2857, NSF 
PHY94-07194 and the CNR. We are grateful to Martin Rees 
and Steinn Sigurdsson for stimulating comments. 
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Fig. 1 . — Upper panels: the fraction of the wind energy radiated by a pulse in the 50-300 keV band as function of pulse duration. 
Parameters: redshift z = 1, r\ m = 30, s e = 0.25, p = 2.5, t w = 20 s, t v = 20 ms (N = 10 3 shells). Only a small fraction of the total 
number of pulses is shown; the density of the points illustrates the pulse duration distribution. The curves shown are log-log space fits 
for the most efficient pulses. The actual values are scattered around the fit. Lower panels: the corresponding PDS (note that f 5 ^Pf 
is plotted in panels Id and 1/), calculated by averaging the PDSs of 100 bursts in each case, all bursts having the same parameters. 
Panels la and lb: optically thin uniform winds, panels lc and Id: optically thin modulated winds; for both cases t]m = 1000 and 
L w = 10 53 ergs _1 . Panels le and If: quasi-optically thick uniform winds (rj M = 300, L w = 10 54 ergs _1 ). Other parameters are 
given in the legend. 




10" 2 10" 1 10° 

f[Hz] 

Fig. 2. — Effect of L w and z on the PDS of optically thin modulated winds. The solid thick line shows the features (i.e. only 
trends, the ordinate values depend on the choice of normalization) of the PDS obtained by BSS98 for real bursts. Other parameters: 
C = 10~ 2 , e e = 0.25, p = 2.5, e B = 10~ 6 , r] m = 30, r] M = 1000, t w = 20 s, t v = 20 ms. Each PDS is the average of 250 power 
spectra of simulated bursts. 
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Fig. 3. — Intensity distribution (panel a), average power density spectrum (panel b) and distribution of Pf around the average PDS 
(panel c) for two models: sin 2 -modulated wind (squares or solid curve) and uniform wind (diamonds or dotted curve). The winds 
have a power-law luminosity distribution and an un-evolving rate density. The mean redshift and dispersion of the bursts brighter 
than 1 7 cm~ 2 s~' are z — 1 .24, cr 2 = 0.91 for the modulated wind model and z — 0.94, a z = 0.84 for the uniform wind model. The 
former model has L m = 7 x 10 52 erg s _1 , tj^j = 1000, Eb = 10~ 6 , ( = 10 -2 , the latter is characterized by L m — 8 x 10 53 erg s~' , 
i]M = 220, Eb = 10 -4 , ( = 1. Other parameters are LM/L m = 100, (3 = 2, r\ m = 30, t w = 20 s, t v = 20 ms, e e = 0.25, p = 2.5, 
for both models. 
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Fig. 4. — Typical 50-300 keV light-curves of the bursts whose PDSs are shown in Figure 3, within the two models for ejection LFs: 
panel a - modulated wind with L w = 10 53 ergs~' , panel b - uniform wind with L w = 10 54 ergs _1 . Both bursts are placed at 
redshift z = 1. Other parameters are as for Figure 3. 



